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ABSTRACT 

Neither  standard  model  SEDs  nor  trancated  standard  model  SEDs  fit  observed  spectra  of  QU  Carinae  with  acceptable 
accuracy  over  the  range  900-3000  A.  Nonstandard  model  SEDs  fit  the  observation  set  accurately.  The  nonstandard 
accretion  disk  models  have  a  hot  region  extending  from  the  white  dwarf  to  R=  1.36f?wD,  a  narrow  intermediate- 
temperature  aimulus,  and  an  isothermal  remainder  to  the  tidal  cutoff  boundary.  The  models  include  a  range  of  M  values 
between  1.0  x  10^’  and  1.0  x  10^^  Mq  yr^^  and  limiting  values  of  A/wq  between  0.6  and  1.2  Mq.  A  solution  with 
TfwD  =1.2  Mq  is  consistent  with  an  empirical  mass-period  relation.  The  set  of  models  agree  on  a  limited  range  of  pos¬ 
sible  isothermal  region  Tiff  values  between  14,000  and  1 8,000  K.  The  model-to-model  residuals  are  so  similar  that  it  is 
not  possible  to  choose  a  best  model.  The  Hipparcos  distance,  6 1 0  pc,  is  representative  of  the  model  results.  The  orbital 
inclination  is  between  40°  and  60°. 

Subject  headings:  novae,  cataclysmic  variables  —  stars:  individual  (constellation  name:  QU  Carinae)  — 
white  dwarfs 


1.  INTRODUCTION 

Cataclysmic  variables  (CVs)  are  semidetached  binary  stars  in 
which  a  late-type  main-sequence  star  loses  mass  onto  a  white  dwarf 
( WD)  by  Roche  lobe  overflow  (Warner  1995).  A  subclass,  the 
NL  systems  (discussed  below)  are  of  interest  because  they  are 
photometrically  stable.  It  is  believed  that  an  analytic  expression 
describes  the  radial  temperature  profile  on  the  accretion  disk, 
leading  to  a  corresponding  spectral  energy  distribution  (SED). 
Tests  using  stellar  SEDs  (Wade  1 988)  produced  unsatisfactory  re¬ 
sults,  but  a  theoretical  model  tailored  to  accretion  disks  (Hubeny 
1990)  is  available  and  might  be  expected  to  lead  to  satisfactory 
results.  In  this  paper  we  show  that  synthetic  SEDs,  built  on  the 
theoretically  prescribed  model,  do  not  fit  the  observed  data,  but 
nonstandard  empirical  models  fit  the  observations  closely. 

In  nonmagnetic  systems  the  mass  transfer  stream  produces  an 
accretion  disk  with  mass  transport  inward  and  angular  momen¬ 
tum  transport  outward  by  viscous  processes.  The  accretion  disk 
may  extend  inward  to  the  WD;  the  outer  boundary  extends  to  a 
tidal  cutoff  limit  imposed  by  the  secondary  star  in  the  steady  state 
case.  If  the  mass  transfer  rate  is  below  a  certain  limit,  the  accretion 
disk  is  unstable  and  undergoes  brightness  cycles  (outbursts),  and 
if  above  the  limit,  the  accretion  disk  is  stable  against  outbursts 
(Shafter  et  al.  1986;  Cannizzo  1993;  Warner  1995;  Osaki  1996). 

’  Based  on  observations  made  with  the  NASA/ESA  Hubble  Space  Telescope, 
obtained  at  the  Space  Telescope  Science  Institute,  which  is  operated  by  the  As¬ 
sociation  of  Universities  for  Research  in  Astronomy,  Inc.,  under  NASA  contract 
NAS5-26555,  and  the  NASA-CNES-CSA  Far  Ultraviolet  Explorer,  which  is 
operated  for  NASA  by  the  Johns  Hopkins  University  under  NASA  contract 
NAS5-32985. 
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The  latter  case  objects  are  called  nova-like  (NL)  systems;  their 
accretion  disks  are  nearly  fully  ionized  to  their  outer  (tidal  cutoff) 
boundary  and  this  condition  suppresses  dwarf  nova  outbursts.  NL 
systems  are  of  special  interest  because  they  are  expected  to  have 
an  accretion  disk  radial  temperature  profile  given  by  an  analytic 
expression  (Frank  etal.  1992,  hereafter  FKR92,  eq.  [5.41]),  thereby 
defining  the  so-called  standard  model. 

Based  on  reports  of  QU  Carinae  as  an  emission-line  object 
(Stephenson  et  al.  1968),  Gilliland  &  Phillips  (1982,  hereafter 
GP82)  undertook  a  spectroscopic  study  to  determine  orbital  char¬ 
acteristics  of  the  system.  They  determined  an  orbital  period  of 
10.9  hr.  No  secondary  star  spectram  was  detected,  indicating  that 
the  spectram  is  dominated  by  light  from  an  accretion  disk  or  the 
primary  star.  Balmer  line  emission  is  relatively  weak.  This  evi¬ 
dence  in  turn  suggested  a  high  rate  of  mass  transfer  and  a  corre¬ 
spondingly  bright  system.  It  also  suggests  that  QU  Car  is  a  NL-type 
system.  We  discuss  known  parameters  of  QU  Car  in  more  detail  in 
§4. 

High-speed  photometry  confirmed  erratic  variations  of  0.1- 
0.2  mag.  on  timescales  of  minutes,  as  reported  by  Schild  (1969). 
Knigge  et  al.  ( 1 994)  used  lUE  spectra  to  study  time  variability  of 
the  UV  resonance  lines.  They  found  that  there  are  short-timescale 
variations  with  a  possible  correlation  with  orbital  period,  but  the 
data  are  insufficient  to  confirm  the  orbital  nature  of  the  variations. 

Hartley  et  al.  (2002)  obtained  Hubble  Space  Telescope  STIS 
spectra  of  QU  Car  at  three  epochs.  Blueshifted  absorption  in  N  v 
and  C IV  provides  evidence  of  an  outflow  with  a  maximum  veloc¬ 
ity  of  «2000  km  s^^  Further  analysis  of  these  spectra  is  in  Drew 
et  al.  (2003).  The  latter  authors  suggest  that  the  QU  Car  secondary 
is  a  carbon  star,  and  that  QU  Car  is  the  most  luminous  NL  system 
known,  possibly  at  a  distance  of  f»2000  pc. 

Our  objective  in  this  study  is  to  use  synthetic  spectra,  based  on 
a  physical  model  of  QU  Carinae,  to  constrain  properties  of  this 
interesting  system. 

2.  INTERSTELLAR  REDDENING 
AND  THE  HIPPARCOS  DISTANCE 

Verbunt  (1987)  determined  E(B-V)  =  Q.\±  0.03  for  QU  Car 
from  the  2200  A  interstellar  medium  (ISM)  feature.  Drew  et  al. 
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(2003)  used  the  STIS  spectra  (the  same  applied  in  this  investiga¬ 
tion)  to  determine  a  neutral  hydrogen  column  density  of  (6  ±  1)  X 
10^0  cm^^.  For  a  standard  gas-to-dust  ratio  of  A(Hi)/Eb-v  = 
5.8  X  10^*  mag^*  cm^^  (Bohlin  et  al.  1978),  the  corresponding 
reddening  is  E(B  —  F)  =  0.1±  0.015,  in  agreement  with  Verbunt 
(1987).  In  §  8  we  model  the  ISM  atomic  and  molecular  hydrogen 
absorption  lines  and  find  a  hydrogen  column  density  of  A(H  i)  = 
2.5x10^®  cm^^,  slightly  lower  than  the  value  obtained  by  Drew 
et  al,  but  within  the  range  of  acceptable  values  for  a  reddening  of 
E{B  —  V  )  =  0.\  (the  Bohlin  et  al.  empirical  law  has  a  very  large 
scatter).  We  adopt  E{B  —  V)  =  O.l. 

The  Elipparcos  parallax  (Perryman  et  al.  1997)  of  1.64  ± 
1.50  (mas),  corresponds  to  a  distance  of  610  pc,  with  1  a  errors 
giving  a  range  between  3 1 8  and  7143  pc.  Thus  the  Elipparcos  par¬ 
allax  does  not  impose  a  strong  constraint. 

3.  THE  FUSE,  STIS,  AND  lUE  SPECTRA 

The  FUSE  and  STIS  spectra  were  obtained  at  a  time  separation 
of  about  3  years.  A  conversion  factor  was  anticipated  to  match  the 
two  spectra,  and  we  found  it  necessary  to  multiply  the  STIS  data 
by  the  factor  1 .263  to  match  the  flux  level  of  the  FUSE  spectrum 
in  the  overlap  region.  Note  also  that  the  STIS  spectrum  used  in 
our  analysis  is  a  combination  of  three  exposures  obtained  over 
an  interval  of  nearly  1  year.  The  data  were  dereddened  to  E{B  — 
F)  =  0.1  (§  2)  after  matching  the  flux  levels. 

The  FUSE  spectrum  of  QU  Car  was  obtained  through  the 
30"  X  30"  LWRS  Large  Square  Aperture,  in  TIME  TAG  mode 
and  was  obtained  on  2003  April  23,  data  set  D 1560 102000  (PI: 
Hartley.  The  data  were  processed  with  CalFUSE  version  3.0.7 
(Dixon  et  al.  2007),  which  automatically  handles  event  bursts. 
Event  bursts  are  short  periods  during  an  exposure  when  high 
count  rates  are  registered  on  one  or  more  detectors.  The  bursts 
exhibit  a  complex  pattern  on  the  detector;  their  cause  is  as  yet 
unknown  (it  has  been  confirmed  that  they  are  not  detector  ef¬ 
fects).  The  main  change  from  previous  versions  of  CalFUSE  is 
that  now  the  data  are  maintained  as  a  photon  list  (the  interme¬ 
diate  data  file  [IDF])  throughout  the  pipeline.  Bad  photons  are 
flagged  but  not  discarded,  so  the  user  can  examine  the  filter  and 
combine  data  without  rerunning  the  pipeline. 

The  spectrum  actually  consists  of  three  exposures,  each  cover¬ 
ing  a  portion  of  a  single  orbit.  These  exposures,  totaling  10,748  s 
of  good  exposiue  time,  were  combined  (using  the  IDF  files).  Con¬ 
sequently,  the  FUSE  spectmm  has  a  diuation  of  about  0.27  of  the 
period  of  QU  Car.  The  spectral  regions  covered  by  the  FUSE  spec¬ 
tral  channels  overlap,  and  these  overlap  regions  are  then  used  to 
renormalize  the  spectra  in  the  SiCl,  LiF2,  and  SiC2  channels  to 
the  flux  in  the  LiFl  channel.  We  combined  the  individual  chan¬ 
nels  to  create  a  time-averaged  spectrum  with  a  linear  0.05  A  dis¬ 
persion,  weighting  the  flux  in  each  output  datum  by  the  exposure 
time  and  sensitivity  of  the  input  exposure  and  channel  of  origin. 
In  this  manner  we  produced  a  final  spectrum  that  covers  the  full 
FUSE  wavelength  range  905-1187  A. 

The  STIS  spectra  used  in  this  study  are  those  obtained  in  2000  by 
Hartley  et  al.  (2002).  The  three  STIS  spectra  were  processed  with 
CALSTIS  version  2.22.  The  HST  STIS  spectra  of  QU  Car  were 
obtained  in  TIME  TAG  operation  mode  using  the  far-ultraviolet 
(FUV)  MAMA  detector.  The  observations  were  made  through 
the  0.2  X  0.2  arcsec^  aperture  with  the  E140M  optical  element. 
Additional  details  are  in  Hartley  et  al.  (2002).  The  STIS  spectra 
have  a  total  duration  of  7500  s,  or  about  0.2  of  the  period  of 
QU  Car.  The  STIS  spectmm,  centered  on  the  wavelength  1425  A, 
consists  of  42  echelle  spectra  that  we  extracted  and  combined  to 
form  a  single  spectmm  extending  irom  about  1160tol710A  and 
with  a  resolution  of  0.1  A. 


lUE  spectra  extend  the  wavelength  coverage  to  3000  A  from 
the  STIS  limit  of  1710  A,  thereby  providing  a  significantly  en¬ 
hanced  constraint  on  a  synthetic  spectmm  model.  The  lUE  archive 
contains  six  SWP  spectra  and  two  LWP  spectra  of  QU  Car  for 
the  dates  1991  June  26-27,  among  others.  Because  of  the  intrin¬ 
sic  variability  of  the  system,  we  have  restricted  our  lUE  data  to 
that  time  interval. 

We  chose  a  single  pair  of  lUE  spectra,  SWP41927  and 
LWP20700,  which  were  taken  in  immediate  sequence  on  1991 
June  27.  The  close  temporal  sequence  minimizes  the  possibility 
of  source  changes  atfecting  the  observations.  We  dereddened 
the  spectra  individually  and  merged  them  with  the  lUEMERGE 
facility,  with  the  “concatenate”  option.  We  used  an  expanded 
ordinate  scale  to  study  the  overlap  of  the  STIS  and  merged  lUE 
spectra  and  found  that  visual  inspection  provided  a  sensitive  test 
of  accurate  superposition.  In  this  way  we  found  that  the  optimum 
divisor  to  superpose  the  merged  spectmm  on  the  scaled  and  de¬ 
reddened  STIS  spectmm  was  0.85. 

Figure  1  shows  the  FUSE  and  the  STIS  spectra  of  QU  Car  be¬ 
fore  dereddening.  Note  the  very  large  number  of  interstellar  lines 
in  the  FUSE  spectmm.  In  Table  1  we  identify  the  prominent  ab¬ 
sorption  lines  of  the  FUSE  spectmm.  The  wavelength  tabulations 
(second  column)  are  the  theoretical  wavelengths.  For  clarity  we 
did  not  include  the  multitude  of  ISM  lines  in  the  table.  A  com¬ 
plete  list  of  ISM  lines  (hydrogen  -I-  metals)  can  be  found  e.g.,  in 
Sembach  (1999)  and  Sembach  et  al.  (2001).  Table  2  identifies  the 
absorption  lines  in  the  STIS  spectmm.  The  ISM  lines  of  the  STIS 
spectmm  have  been  analyzed  in  Drew  et  al.  (2003),  and  the  broad 
absorption  lines  (C  m  1 1 1 76,  N  v  2 1 240, 0  v  1 1 37 1 ,  Si  iv  2 1 398, 
C IV  21549,  He  ii  21640)  have  been  studied  in  detail  for  each  STIS 
exposure  separately  in  Hartley  et  al.  (2002).  The  broad  lines  are 
seen  in  emission  with  a  velocity  broadening  of  several  thousand 
km  s^^;  an  absorption  component  is  superposed  with  a  velocity 
broadening  of  a  few  hundred  km  s^^ 

4.  PARAMETERS  ADOPTED  INITIALLY 

The  orbital  period  of  10.9  hr  was  determined  by  GP82,  using 
radial  velocities  from  the  224630-4660  emission  feature  and  the 
24686  He  ii  emission  line.  There  is  no  detectable  orbital  light  mod¬ 
ulation  (Schild  1969;  GP82),  suggesting  an  inclination  i  <  60°. 
The  presence  of  doubled  emission  lines  (GP82)  suggests  i>  30°. 

Based  on  the  main-sequence  mass-radius  relation  of  Lacy 
(1977)  and  the  usual  period-mean  density  relation  (Robinson 
1976),  GP82  obtained  an  approximate  M2(=Ms)  =  1.3  Mq. 
Note  that  this  determination  does  not  depend  on  an  explicit 
knowledge  of  i.  GP82  spectrophotometry  sets  a  limit  of  «1% 
on  the  depth  of  any  absorption  feature  from  the  secondary  star  in 
the  QU  Car  optical  spectmm.  Comparison  with  the  « 1 0%  depth 
of  absorption  features  in  G-K  standards  (appropriate  to  the  sec¬ 
ondary  star)  establishes  that  the  disk  and  primary  of  QU  Car  must 
be  radiating  at  ^10  Lq.  The  strong  He  ii  24686  emission,  weaker 
He  II 24542  emission,  and  the  absence  of  He  i  lines  at  243 87, 2447 1 , 
and  24713  suggests  a  temperature  in  the  line-emitting  region  in 
excess  of  20,000  K.  This  is  consistent  with  a  fully  ionized  accre¬ 
tion  disk  expected  for  a  NL  system.  GP82  determine  a  minimum 
distance  to  QU  Car  of  500  pc,  based  on  mb  =  11.2  for  the  sys¬ 
tem  and  standard  assumptions  for  interstellar  absorption.  From 
the  My  —  Porb  relation  of  Warner  (1995,  Fig.  4.16),  with  My  = 
2.4  and  interstellar  absorption  of  0.3  mag.  (§  2),  we  obtain  a 
minimum  distance  that  agrees  with  the  GP82  value. 

A  coarse  estimate  of  the  mass  transfer  rate  is  possible  using  the 
study  by  Smak(  1989).  From  Smak’s  Figure  2  we  obtain  M  ^  5  x 
10^’  Mq  yr^*.  Later,  Smak  (1994)  slightly  revised  the  calculated 
mass  transfer  rates.  This  value  is  roughly  consistent  with  the  mass 
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Fig.  1 . — FUSE  plus  STIS  combined  spectra  with  line  identifications.  The  STIS  spectrum  has  been  multiplied  by  1 .263  to  match  the  flux  level  of  the  FUSE  spectrum  in 
the  overlap  region. 


transfer  rate-orbital  period  relation  of  Patterson  (1984,  Fig.  21). 
Drew  et  al.  (2003)  find  that  the  QU  Car  C/Fle  abundance  may  be 
as  high  as  0.06,  an  order  of  magnitude  higher  than  the  solar  ratio, 
and  the  C/0  abundance  is  estimated  to  be  greater  than  1 .  Based  on 
a  power- law  fit  to  merged /C/E  spectra,  and  a  possible  cutoff  of  the 
spectmm  at  228  A  (the  He  ii  series  limit).  Drew  et  al.  (2003)  sug¬ 
gest  a  mass  transfer  rate  of  M  =  8 . 0  x  lO-*M0yr-‘if  QU  Car  is 
at  a  distance  of  500  pc  and  M  =  1 .0  x  10^®  Mq  yr^*  if  the  dis¬ 
tance  is  2000  pc. 


GP82  cite  the  uncertainty  in  the  value  and  the  lack  of 
knowledge  of  i,  and  so  do  not  derive  a  value  of  Mwd-  However, 
some  constraint  on  A/wd  is  possible  as  shown  by  the  following 
discussion. 

From  the  standard  expressions  for  radial  velocity  amplitude 
(Smart  1947  p.  360,  eq.  60),  we  obtain 


2tt  sin  i 
A/wd-P 


(1) 
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TABLE  1 
FUSE  Lines 


Wavelength 

(A) 

Line  Identification 

Origin 

FWHM 

(A) 

Comments 

See  text . 

H  I 

C,  ISM 

See  text . 

H  II 

C,  ISM 

See  text . 

D  I 

C,  ISM 

See  text . 

0  I 

C,  ISM 

923.0 . 

N  IV  921.46-924.91 

s 

4.0 

Unresolved  multiplet 

933.3 . 

S  VI  933.38 

s 

3.0 

Affected  by  H2 

944.6 . 

S  VI  944.70 

s 

2.5 

Affected  by  H2 

950.5 . 

P  IV  950.66 

s? 

Heavily  affected  by  H2 

951.05 . 

N  I  951.12 

C,  ISM 

0.1 

953.35 . 

N  I  953.42 

C,  ISM 

0.1 

953.60 . 

N  I  953.66 

C,  ISM 

0.1 

953.97 . 

N  I  953.97  +  954.10 

C,  ISM 

0.3 

Unresolved 

961.00 . 

P  II  961.04 

C,  ISM 

0.1 

962.10 . 

P  II  962.12 

C,  ISM 

0.1 

963.75 . 

P  II  963.80 

C,  ISM 

0.1 

Affected  by  H2 

963.95 . 

N  I  963.99 

C,  ISM 

0.1 

964.57 . 

N  I  964.63 

C,  ISM 

0.1 

977.00 . 

C  III  977.03 

C,  ISM 

0.15 

977.00 . 

C  III  977.03 

S? 

2.0 

Affected  by  H2 

989.80 . 

N  III  989.77  +  Si  ii  989.89 

C,  ISM 

Affected  by  H2 

1031.80 . 

0  VI  1031.65 

s 

3.5 

1036.34 . 

C  II  1036.34  doublet 

ISM 

Affected  by  H2 

1037.2? . 

0  VI  1037.30 

S 

3.0 

Heavily  affected  by  H2 

1048.20 . 

Ar  I  1048.22 

ISM 

0.1 

1055.25 . 

Fe  II  1055.27 

ISM 

0.1 

1062.6? . 

S  IV  1062.66 

S 

Heavily  affected  by  H2 

1063.15 . 

Fe  II  1063.18 

ISM 

0.1 

1066.65 . 

Ar  I  1066.66 

ISM 

0.1 

1073.2 . 

S  IV  1073.1  +  1073.5 

S 

2.0 

Unresolved  doublet 

1081.85 . 

Fe  II  1081.87 

C,  ISM 

Affected  by  H2 

1083.97 . 

N  II  1083.99 

C,  ISM 

0.15 

1096.85 . 

Fe  II  1096.88 

C,  ISM 

0.1 

Affected  by  H2 

1112.00 . 

Fe  II  1112.03 

C,  ISM 

0.1 

1118.05 . 

S  IV,  VI  1117.9  +  P  V  1117.98 

s 

1.7 

1121.97 . 

Fe  II  1121.97 

C,  ISM 

0.1 

1122.40 . 

Si  IV  1122.48 

s 

1.5 

1128.20 . 

P  V  1128.01+  Si  IV  1128.33 

s 

1.5 

1133.70 . 

Fe  II  1133.67 

C,  ISM 

0.1 

1134.15 . 

N  I  1134.17 

A,  C,  ISM 

0.1 

1134.45 . 

N  I  1134.42 

A,  C,  ISM 

0.1 

1135.00 . 

N  I  1134.98 

A,  C,  ISM 

0.15 

1142.40 . 

Fe  II  1142.37 

C,  ISM 

0.1 

1143.25 . 

Fe  II  1143.23 

C,  ISM 

0.1 

1144.95 . 

Fe  II  1144.94 

C,  ISM 

0.1 

1152.85 . 

P  II  1152.82 

C,  ISM 

0.1 

1168.75 . 

He  I  1168.9 

A?,  S? 

0.3 

Emission 

1172.50 . 

S  I  1172.55? 

C,  ISM 

0.3 

1173.00 . 

0  III  1173.37 

S? 

0.3 

1175.95 . 

C  III  1174.9-1176.4 

S 

1.7 

Shifted  by  +0.35  A 

Note.  — A  =  possibly  including  air  glow  contamination;  C  =  circumbinary  material;  ISM  =  ISM;  S  =  source. 


where  we  assume  zero  orbital  eccentricity;  the  mass  ratio  q  — 
A/j/A/wd,  i  is  the  orbital  inclination,  is  the  WD  radial  ve¬ 
locity  amplitude,  P  is  the  orbital  period,  and  D  is  the  component 
separation.  From  Kepler’s  third  law. 


47r^  1 


(2) 


From  these  two  equations  we  have 


=  GMs 


P  sin  i 

2tt  Kwd 


(3) 


GP82  determine  a  value  Kwd  =  115  km  s^  *  from  He  ii  emission 
lines  (insecure  because  emission  lines  may  not  trace  the  WD  mo¬ 
tion  accurately).  GP82  derive  an  orbital  period  of  0.4542  days. 
With  these  parameters  determined  by  observation,  an  adopted 
value  of  i  and  an  assigned  value  of  determines  D  from  equa¬ 
tion  (3),  q  from  equation  (1),  and  A/wp  from  Afwo  =  MJq.  Thus, 
a  given  i  has  an  associated  contour  in  the  mass-mass  plane  that 
determines  AT^d  values  corresponding  to  assigned  values. 

Figure  2  illustrates  the  principle.  The  WD  axis  is  bounded 
above  by  the  Chandrasekhar  mass.  The  diagonal  line  represents  the 
case  of  equal  WD  and  secondary  star  masses;  systems  represented 


TABLE  2 
STIS  Lines 


Wavelength  FWHM 

(A)  Line  Identification  Origin  (A)  Comments 

1168.59 .  N I  1168.54  C,  ISM  0.1 

1171.08 .  N I  1171.08  C,  ISM  0.15 

1175.7  .  Cm  S  Unresolved  multiplet 

1188.83  .  C  I  1188.83  C,  ISM  0.1 

1190.23  .  Sin  1190.20  C,  ISM  0.07 

1190.29  .  C  I  1190.25  C,  ISM  Unresolved 

1190.43  .  Sin  1190.42  C,  ISM  0.13 

1193.05 .  C  I  1193.03  C,  ISM  0.07 

1193.30  .  Sin  1193.29  C,  ISM  0.16 

1197.21 .  Mn  II  1197.18  C,  ISM  0.05 

1199.43  .  Mn  II  1199.39  C,  ISM  0.05 

1199.56  .  N I  1199.55  C,  ISM  0.16 

1200.23  .  N I  1200.22  C,  ISM  0.14 

1200.72  .  N I  1200.71  C,  ISM  0.13 

1201.15 .  Mnn  1201.12  C,  ISM  0.05 

1206.53  .  Si  in  1206.50  C,  ISM  0.12 

1215.7  .  Hi  1215.66  S,  C,  ISM  12 

1237.09 .  Ge  n  1237.06  C,  ISM  0.06 

1238.94  .  N  V  1238.82  S  2.2 

1239.96 .  Mg  n  1239.93  C,  ISM  0.05 

1240.42  .  Mg  n  1240.39  C,  ISM  0.05 

1243.00 .  N  V  1242.80  S  1.7 

1250.62 .  Sn  1250.58  C,  ISM  0.11 

1253.83  .  Sn  1253.81  C,  ISM  0.12 

1259.53  .  Sn  1259.52  C,  ISM  0.12 

1260.45 .  Si  II  1260.42?  ?  0.26  Asymmetric 

1260.77 .  C  I  1260.74  C,  ISM  0.06 

1266.4 .  C  I  1266.42  S  1.6  Shallow 

1272.7  .  0  n  1272.23  S  1.5  Shallow 

1277.29 .  Znn  1277.31?  C,  ISM  0.08 

C  I  1277.25  +  1277.28?  C,  ISM  0.08 

1277.57? .  C  I  1277.51  +  1277.55?  C,  ISM  0.12  Very  shallow,  unresolved 

1280.18  .  C  I  1280.14  C,  ISM  0.05  Shallow 

1301.91 .  Pn  1301.87  C,  ISM  0.05 

1302.19  .  0  I  1302.17  C,  ISM  0.20 

1304.39 .  Si  II  1304.37  C,  ISM  0.14 

1304.85  .  0  I  1304.86  C,  ISM  0.20  Very  shallow 

1306.02 .  0  I  1306.03  C,  ISM  0.20  Very  shallow 

1317.25 .  Ni  II  1317.22  C,  ISM  0.09  Shallow 

1328.87 .  C  I  1328.83  C,  ISM  0.06  Shallow 

1329.14 .  C  I  1329.12  C,  ISM  0.06  Very  shallow 

1334.56  .  C  n  1334.53  C,  ISM  0.21 

1335.73  .  Cn  1335.71  C,  ISM  0.14 

1347.28 .  Cl  I  1347.24  C,  ISM  0.05  Shallow 

1355.65 .  0  I  1355.60  C,  ISM  0.05  Very  shallow 

1358.81  .  Cun  1358.77  C,  ISM  0.05  Very  shallow 

1370.17  .  Ni  II  1370.14  C,  ISM  0.07  Shallow 

1371.6  .  Ov  1371.30  S  1.6 

1393.95  .  Si  IV  1393.76  S  1.1 

1403.95  .  Si  IV  1402.77  S  1.2 

1420.7  .  ?  S?  1.7  Seen  only  in  one  exposure 

1414.43  .  Ga  ii  1414.40  C,  ISM  0.05  Very  shallow 

1454.89 .  Ni  n  1454.84  C,  ISM  0.04  Very  shallow 

1526.73  .  Si  II  1526.71  C,  ISM  0.17 

1532.58 .  P  n  1532.53  C,  ISM  0.05  Very  shallow 

1548.18  .  Civ  1548.19  S  1.0 

1550.85  .  C  IV  1550.77  S  1.0 

1560.36 .  C  I  1560.31  C,  ISM  0.06 

1560.74  .  C  I  1560.71  C,  ISM  0.07  Very  shallow 

1640.52 .  He  ii  1640.47  S  1.0 

1656.98 .  C  I  1656.93  C,  ISM  0.07 

1657.43  .  C  I  1657.38  C,  ISM  0.04  Very  shallow 

1657.96  .  C  I  1657.91  C,  ISM  0.04  Very  shallow 

1670.81  .  All  1670.79  C,  ISM  0.13 


QU  CAR 
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WD  mass{M0) 

Fig.  2. — Mass-mass  plot  showing  regions  where  a  system  model  is  permitted. 
The  region  above  the  diagonal  line  is  excluded  by  dynamical  instability.  The  region 
below  i  —  60°  is  excluded  by  the  absence  of  observed  eclipses.  See  the  text  for 
details. 


by  points  above  this  line  are  excluded  by  dynamical  instability. 
( Mass  transfer  from  a  more  massive  secondary  causes  the  period 
to  shorten,  the  separation  of  components  and  the  secondary  com¬ 
ponent  Roche  lobe  to  shrink,  and  so  accelerates  the  mass  transfer.) 
The  value  of  i  cannot  be  greater  than  about  60°  to  avoid  eclipses. 
No  models  with  i  as  small  as  30°  are  permitted  (for  <7=1.0  the 
masses  of  both  stellar  components  then  are  larger  than  the 
Chandrasekhar  mass  and  for  ^<1.0  the  WD  mass  becomes 
even  larger).  A  A4  >  1.30  A/q  is  not  permitted  because  the  sec¬ 
ondary  must  lie  on  or  above  the  ZAMS  (zero-age  main-sequence) 
mass-radius  relation  (Lacy  1977).  The  empirical  mass-period  re¬ 
lation  of  Warner  (1995)  eq.  [2.100]  gives  M,  =  1.29  Mq.  A  value 
Mg  =  \  .20  Mq,  with  A/wd  =  1  -20  Mq  is  permitted  with  /  =  38° . 
The  uncertainty  in  leads  us  to  adopt  i  =  40°,  Afwo  =  Ms  = 
1 .20  Mq  as  an  appropriate  case  for  simulation. 

5.  CALCULATION  OF  SYSTEM  MODELS 

Our  calculation  of  system  models  uses  the  BINSYN  suite 
(Linnell  <&  Hubeny  1996;  Linnell  et  al.  2007a).  The  word  “model,” 
as  used  in  this  section,  has  two  context-dependent  meanings.  In 
one  context  it  refers  to  a  listing  of  the  run  of  physical  parameters 
in  an  annulus,  or  to  properties  of  the  collection  of  annuli  compris¬ 
ing  an  entire  accretion  disk.  In  a  different  context  it  refers  to  the 
array  of  output  synthetic  spectra  from  BINSYN.  In  the  latter  con¬ 
text,  orbital  phase-dependent  synthetic  spectra  separately  of  the 
system,  the  WD,  the  secondary  star,  the  accretion  disk,  and  the  ac¬ 
cretion  disk  rim,  all  produced  in  a  single  run,  constitute  a  BIN S  YN 
model  of  a  CV  system.  The  synthetic  spectra  include  effects  of 
eclipse  and  irradiation. 

A  BINSYN  accretion  disk  model  uses  an  array  of  annulus  mod¬ 
els  calculated  with  the  program  TLUSTY  (Hubeny  1988;  Hubeny 
1990;  Hubeny  <fe  Lanz  1995;  Hubeny  <&  Hubeny  1998)  and  cor¬ 
responding  synthetic  spectra  calculated  with  program  SYNSPEC 
( Hubeny  et  al.  1 994).  Table  3  is  an  illustration  of  an  accretion  disk 
model  for  A/wd  =  L2  Mq  and  M  =  6.0  x  10^^  Mq  yr^^  The 
26  TLUSTY  annulus  models  cover  the  interval  from  the  WD  equa¬ 
tor  to  the  tidal  cutoff  boundary.  The  cutoff  boundary  is  calculated 
from  Warner  (1995,  p.  57,  eq.  [2.61]).  Note  that  the  annulus  radii 
in  the  first  column  are  measured  in  units  of  a  zero-temperature 
WD.  Successive  columns  provide  some  of  the  data  available  in  the 
TLUSTY  output  for  a  given  annulus  model;  each  row  describes 
a  separate  annulus  model.  The  annuli  all  are  H-He  models  (i.e.. 


TABLE  3 

Properties  of  Accretion  Disk  with  Mass  Transfer  Rate 
M  —  6.0  X  10“^  Mq  yr“^  and  WD  Mass  of  1.20  Mq 


r/rwD.o 

TcS 

mo 

T 

logs 

^Ross 

1.36 . 

334,084 

2.213E5 

4.534E6 

7.88 

6.92E7 

8.73E4 

2.00 . 

299,540 

2.701E5 

4.22 1E6 

7.62 

1.21E8 

1.02E5 

3.00 . 

242,217 

2.640E5 

3.386E6 

7.31 

1.99E8 

9.93E4 

4.00 . 

203,586 

2.458E5 

2.802E6 

7.08 

2.76E8 

9.34E4 

5.00 . 

176,588 

2.284E5 

2.394E6 

6.89 

3.56E8 

8.80E4 

6.00 . 

156,665 

2.133E5 

2.097E6 

6.75 

4.38E8 

8.36E4 

7.00 . 

141,312 

2.004E5 

1.871E6 

6.62 

5.20E8 

7.99E4 

8.00 . 

129,082 

1.894E5 

1.693E6 

6.51 

6.04E8 

7.68E4 

9.00 . 

119,082 

1.798E5 

1.550E6 

6.42 

6.90E8 

7.43E4 

10.00 . 

110,733 

1.714E5 

1.432E6 

6.33 

7.77E8 

7.22E4 

12.00 . 

97,540 

1.572E5 

1.247E6 

6.18 

9.54E8 

6.88E4 

14.00 . 

87,537 

1.460E5 

1.110E6 

6.06 

1.14E9 

6.63E4 

16.00 . 

79,657 

1.366E5 

1.004E6 

5.95 

1.32E9 

6.45E4 

18.00 . 

73,268 

1.288E5 

9.197E5 

5.85 

1.51E9 

6.3 1E4 

20.00 . 

67,967 

1.221E5 

8.505E5 

5.77 

1.71E9 

6.21E4 

24.00 . 

59,649 

1.110E5 

7.439E5 

5.62 

2.11E9 

6.08E4 

30.00 . 

53,387 

1.024E5 

6.334E5 

5.50 

2.73E9 

6.02E4 

35.00 . 

45,436 

9.088E4 

5.682E5 

5.32 

3.27E9 

6.03E4 

45.00 . 

37,854 

7.930E4 

4.782E5 

5.12 

4.37E9 

6.20E4 

60.00 . 

30,685 

6.767E4 

3.957E5 

4.89 

6.11E9 

6.61E4 

80.00 . 

24,852 

5.762E4 

3.304E5 

4.66 

8.53E9 

7.35E4 

100.00 . 

21,092 

5.08 1E4 

2.891E5 

4.49 

I.IOEIO 

8.21E4 

120.00 . 

18,440 

4.58 1E4 

2.601E5 

4.34 

1.37E10 

9.15E4 

140.00 . 

16,458 

4.195E4 

2.383E5 

4.22 

1.64E10 

1.02E5 

160.00 . 

14,911 

3.885E4 

2.211E5 

4.11 

1.91E10 

1.12E5 

180.00 . 

13,667 

3.63 1E4 

2.069E5 

4.01 

2.19E10 

1.22E5 

Note. — Each  line  in  the  table  represents  a  separate  annulus.  A  Shakura  & 
Sunyaev  (1973)  viscosity  parameter  a  =  0.1  was  used  in  calculating  all  annuli. 
'■wD,o  is  the  radius  of  a  zero-temperature,  carbon,  Hamada-Salpeter  WD.  See  the 
text  for  additional  comments. 

these  are  the  only  explicit  atoms;  the  next  28  elements  are  treated 
implicitly;  see  the  TLUSTY  manual  for  details*).  Convection  was 
ignored.  The  models  converged  for  all  of  the  Table  3  annuli.  In 
some  instances,  for  other  A7  values,  the  outermost  annuli  failed  to 
converge  and  we  used  so-called  gray  models  (see  the  TLUSTY 
manual). 

The  niQ  column  lists  the  mass  per  unit  area  (g  cm^^)  between 
the  central  plane  and  the  upper  boundary.  Calculation  of  an  an¬ 
nulus  model  specifies  an  initial  uppermost  layer  in  g  cm^^.  Our 
calculation  used  10^^  g  cm^^  for  the  initial  layer.  An  annulus 
model  typically  may  have  70  depth  points.  7],  is  the  central  plane 
temperature,  log  g  is  the  value  of  that  parameter  at  optical  depth 
near  0.7,  Zfj  is  the  height,  in  cm,  of  the  optical  depth  0.7  level 
above  the  central  plane,  and  is  the  Rosseland  optical  depth  at 
the  central  plane. 

SYNSPEC  is  used  to  calculate  a  synthetic  spectrum  for  each 
TLUSTY  annulus,  with  specified  spectral  resolution.  The  line 
list  includes  data  for  the  first  30  periodic  table  elements.  Next, 
BINSYN  is  used  with  its  specified  number  of  annuli  (different 
from  the  set  of  TLUSTY  annuli)  to  calculate  the  various  output 
synthetic  spectra  by  integration  over  the  array  of  BINSYN  annuli 
for  the  accretion  disk  and  by  corresponding  integrations  over  the 
stellar  photospheric  grids  for  the  latter  objects.  Synthetic  spectra 
for  the  individual  BINSYN  annuli  follow  by  interpolation  among 
the  array  of  SYNSPEC  spectra  of  the  TLUSTY  annuli.  It  is  at 
this  level  that  information  is  incorporated  on  orbital  inclination, 
eclipse  effects,  radial  velocity  effects,  etc.,  as  needed  for  individual 
photospheric  segments  of  all  objects  in  the  system. 

*  See  http;//nova.astro.umd.edu. 
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Calculation  of  a  synthetic  spectrum  corresponding  to  a  given 
M  model  is  a  computationally  intensive  process  and  compromises 
are  necessary  in  deciding  how  many  TLUSTY  annulus  models  to 
calculate  to  achieve  acceptable  accuracy  in  the  final  system  syn¬ 
thetic  spectmm.  We  believe  our  representational  accuracy  is  high, 
but  we  cannot  be  sure  that  the  fit  in,  e.g.,  Figure  8  could  not  be 
improved  with  a  finer  grid  of  annulus  synthetic  spectra  for  the 
TLUSTY  annulus  models  and  a  larger  number  of  BINSYN  annuli. 
We  note  that  this  entire  process  must  be  repeated  for  each  value  of 
M  and  in  turn  for  each  value  of  A/wd- 

In  general,  the  WD  synthetic  spectmm  contribution  to  the  sys¬ 
tem  synthetic  spectmm  is  a  strong  function  of  the  adopted  WD 
Tts-  Since  the  WD  radius  varies  with  the  adopted  WD  TeS,  it 
makes  a  difference  whether  we  use  the  zero-temperature  radius 
or  the  corrected  radius  in  calculating  annulus  radii,  in  BINSYN, 
measured  in  units  of  the  WD  radius.  We  must  use  the  corrected 
radius  with  BINSYN,  designated  by  the  symbol  since  the 
output  spectra  must  describe  the  WD  contribution  for  its  ac¬ 
tually  adopted. 

The  WD  Tgff  likely  will  be  unknown  in  advance  of  the  model 
calculation.  We  determine  what  annulus  models  to  calculate  with 
TLUSTY  by  using  a  fixed  set  of  multiples  of  the  zero-temperature 
WD  radius;  we  stress  that,  if  we  were  to  use  the  fixed  multiples 
procedure  with  a  variable  WD  radius,  it  would  be  nearly  prohib¬ 
itive  to  recalculate  the  entire  array  of  TLUSTY  annuli  with  each 
new  approximation  to  an  adopted  WD  Such  iterative  recal¬ 
culation  is  unnecessary;  BESfSYN  temperature-wise  interpolation 
among  the  fixed  (with  rwo  =  ^wd,o)  TLUSTY  annulus  models 
leaves  the  accretion  disk  synthetic  spectmm  essentially  unaffected 
as  the  WD  Teg-dependent  radius  changes.  We  adopt  solar  com¬ 
position  for  the  annulus  models.  In  the  present  application  the 
WD  contribution  is  too  small  to  matter,  but  it  is  important  to  retain 
the  distinction,  as  this  paper  does. 

Our  initial  model,  at  all  mass  transfer  rates,  adopted  a  standard 
model  (FKR92;  Linnell  et  al.  2007a).  A  standard  model  assumes 
the  accretion  disk  is  in  hydrostatic  equilibrium,  that  the  accretion 
disk  material  is  very  nearly  in  Keplerian  rotation,  and  that  a  sin¬ 
gle  viscosity  parameter  applies  to  the  entire  accretion  disk.  All 
of  our  TLUSTY  annulus  models  adopted  a  viscosity  parameter 
(Shakura  &  Sunyaev  1973)  a  =  0.1. 

6.  QU  CARINAE  MODELS 
WITH  Mwd  =  1.20  Mq,  9  =  1.00 

The  long  orbital  period  and  the  weak  Balmer  line  emission, 
suggesting  a  NL-type  system  (GP82),  in  turn  indicate  a  hot  WD. 
The  radius  of  a  1 .2  Mq  WD  is  slightly  sensitive  to  the  photospheric 
temperature.  Table  12C-H-S.dat  of  Panei  et  al.  (2000)  for  homo¬ 
geneous  Hamada-Salpeter  models,  lists  a  radius  rwo.o  =  5.6  x 
10^^  Rq  for  a  zero-temperature,  carbon,  1 .2  Mq  WD.  Based  on  an 
adopted  55,000  K  WD,  we  used  Figure  4a  of  Panei  et  al.  (2000)  to 
determine  a  corrected  radius  of  6.7  x  10  ^  Rq.  We  will  find  that 
the  WD  makes  an  insignificant  contribution  to  the  system  lumi¬ 
nosity;  the  only  impact  the  WD  Egg  has  is  indirect — temperatures 
in  the  accretion  disk  at  specified  multiples  of  the  WD  radius  vary 
with  the  Tgg-dependent  WD  radius. 

A  smaller  assumed  WD  Tgg  produces  a  smaller  and  a  cor¬ 
respondingly  larger  accretion  disk  T^g  in  the  innermost  rings 
( FKR92)  since  they  then  are  in  a  deeper  potential  well.  However, 
this  effect  is  smaller  than  the  temperature  effect,  on  all  annuli,  of 
different  assumed  mass  transfer  rates.  For  uniformity  in  compar¬ 
ing  different  assumed  mass  transfer  rates  we  arbitrarily  adopt  a 
fixed  WD  Teg  =  5 5 , 000  K.  The  tidal  cutoff  radius  of  the  accretion 
disk  is  1  50Rwd  for  the  adopted  WD  Tgg,  and  also  equals  0.998  ^0. 
(Contrast  with  the  same  tidal  cutoff  radius,  measured  in  units 


Fig.  3. — Standard  model  fit  to  FUSE  +  STIS  +  lUE  spectra.  The  synthetic 
spectrum  is  a  continuum  spectrum  +  H  and  He  n  lines.  The  mass  transfer  rate  is  M  — 
6.0  X  10“^  Mq  yr“\  and  the  WD  mass  is  1.20  Mq.  The  synthetic  spectrum  spectral 
gradient  is  too  large.  Note  that  the  ordinate  scale  is  logarithmic. 

of  the  zero-temperature  WD  radius,  in  the  last  line  of  Table  3). 
The  BINSYN  model  uses  45  division  circles  on  the  accretion 
disk,  with  the  first  division  at  the  WD  equator,  the  next  two  at 
1.1806Rwd,  and  1.361  1Rwd>  and  the  last  at  the  tidal  cutoff 
boundary. 

We  tested  models  with  M  =  0.1,  0.3,  0.6,  1,  3,  6,  and  10  x 
10^^  Mq  yr^*  In  the  first  case  a  standard  model  produces  an  ac¬ 
cretion  disk  temperature  profile  whose  outer  annuli  are  cooler 
than  6000  K.  This  model  is  unstable  (Smak  1982)  and  must  be 
rejected.  In  all  remaining  cases,  the  standard  model  produces  an 
accretion  disk  profile  with  a  spectral  gradient  that  is  too  large. 

Figure  3  illustrates  the  problem  for  M  =  6.0x  lO-’M0yr-‘; 
other  values  of  M  produce  comparable  plots  after  applying  ap¬ 
propriate  normalizing  factors.  In  all  cases  in  this  paper,  normal¬ 
izing  factors  have  been  determined  by  eye  estimates.  In  Figure  3 
the  peak  fiux  of  the  WD,  near  1 000  A,  is  less  than  0.0 1  of  the  sys¬ 
tem  flux;  the  WD  contribution  to  the  system  spectmm  is  negli¬ 
gible.  Note  that  the  synthetic  spectmm  consists  of  a  continuum 
spectmm  plus  absorption  lines  of  H  and  He  ii.  (He  ii  but  not  He  i 
because  He  ii  is  a  single  electron  ion.)  The  spectral  gradient  prob¬ 
lem  persists  if  i  is  arbitrarily  changed  to  50°.  For  reasons  discussed 
in  §  10,  the  M  =  1 .0  x  IO^^Mq  yr^*  rate  is  a  maximum,  so  our 
range  of  M  values  cover  the  full  range  of  possible  choices  for 
this  M^y)' 

We  conclude  that,  for  this  WD  and  secondary  star  mass,  no  stan¬ 
dard  model  SED  provides  an  acceptable  fit  to  the  observational 
data.  We  believe  the  spacing  between  our  M  values  is  small  enough 
to  extrapolate  our  conclusion  to  any  M  value  not  actually  calculated. 

Tmncated  standard  models  are  of  interest  because  the  inner 
accretion  disk  can  be  elided  by  radiation  from  a  hot  WD  or  by 
interaction  with  a  magnetic  field  associated  with  the  WD.  We  ig¬ 
nore  whether  either  process  is  likely  for  QU  Car  and  restrict  con¬ 
sideration  to  whether  a  tmncated  accretion  disk  can  provide  a  SED 
that  fits  the  observations  with  acceptable  accuracy. 

A  tmncated  accretion  disk  model,  with  the  residual  disk  on  the 
standard  model  temperature  profile,  should  have  a  lower  spectral 
gradient  and  so  is  a  possible  alternative.  We  have  calculated  tmn¬ 
cated  models  for  the  last  six  mass  transfer  rates  listed  above  (the 
model  with  the  lowest  rate  of  mass  transfer  failed  because  of  its 
temperature  behavior  at  large  radii).  Tmncation  cuts  off  the  annuli 
that  contribute  most  to  the  far  UV.  In  the  six  largest  M  cases,  the 
best  tmncation  compromise  either  still  had  a  too  large  spectral  gra¬ 
dient  or  failed  to  represent  the  FUSE  spectral  region  with  accept¬ 
able  fidelity. 
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Fig.  4. — Truncated  model  fit  to  observed  spectra.  The  mass  transfer  rate  is  M  = 
6.0  X  10“^  Mq  yr~\  the  WD  mass  is  1.20  Mq,  and  the  accretion  disk  has  been 
truncated  at  /?tn,n  =  22.5i?wD-  The  spectral  gradient  is  a  better  fit  than  Fig.  3  but 
still  unacceptable  and  the  synthetic  spectrum  fails  to  fit  the  FUV  accurately. 

Figure  4  illustrates  the  fit  for  A/  =  6.0  x  lO^^Afg  yr^',  and  a 
truncation  radius  of  Rtmnc  =  22.5Rwd-  The  overall  fit  is  an  im¬ 
provement  over  Figure  3,  but  the  fit  to  the  FUSE  spectrum  is  un¬ 
satisfactory  and  the  spectral  gradient  still  is  too  large.  A  larger 
truncation  radius  would  improve  the  spectral  gradient  problem 
but  make  the  fit  to  the  FUSE  spectrum  worse.  A  smaller  trunca¬ 
tion  radius  has  a  still  more  discrepant  spectral  gradient.  Other  val¬ 
ues  of  M  are  unsatisfactory  for  similar  reasons. 

Having  demonstrated  that  no  standard  model  SED  or  tmncated 
standard  model  SED  fits  the  observational  data,  we  consider  non¬ 
standard  model  SEDs.  The  too-large  spectral  gradient  can  be  re¬ 
duced  with  a  cooler  (average)  accretion  disk.  However,  the  flux 
maximum  in  the  FUSE  spectrum  near  1000  A  requires  the  pres¬ 
ence  of  a  hot  region  on  the  accretion  disk.  This  suggests  reducing 
the  temperature  of  the  disk  model  at  intermediate  radii. 

It  is  important  to  note  that  a  modification  of  the  temperature 
profile  removes  the  direct  connection  to  a  specific  mass  transfer 
rate.  Not  only  is  the  empirical  profile  disconnected  from  the  stan¬ 
dard  model,  the  total  flux  from  the  empirical  model  differs  from 
the  flux  of  the  parent  standard  model.  We  conducted  very  exten¬ 
sive  preliminary  experiments  and  found  that  a  steep  drop  in  the 
accretion  disk  temperature,  followed  by  a  much  lower  tempera¬ 
ture  gradient  than  in  the  standard  model  gave  a  fairly  good  fit  to 
the  observed  spectra.  However,  we  determined  that  further  im¬ 
provement  was  possible. 

It  is  desirable  to  preserve  model  connection  to  some  mass  trans¬ 
fer  rate.  On  the  standard  model,  half  of  the  potential  energy  liber¬ 
ated  by  mass  falling  from  the  El  point  (located  on  the  secondary 
star  photosphere)  to  the  surface  of  the  WD  is  converted  to  emitted 
radiation  by  the  accretion  disk  (FKR92).  The  remaining  half 
either  is  liberated  in  a  boundary  layer,  goes  into  spinning  up  the 
WD,  or  is  deposited  elsewhere,  as  in  a  wind.  Based  on  energy  con¬ 
siderations,  a  nonstandard  model  accretion  disk  still  can  be  asso¬ 
ciated  with  a  mass  transfer  rate  if  roughly  half  of  the  total  liberated 
potential  energy  appears  as  flux  from  the  accretion  disk. 

With  that  constraint  in  mind,  we  performed  further  tests.  In  the 
standard  model  the  maximum  accretion  disk  temperature  occurs 
at  R(Tn,ax)  =  1.3611Rwd-  The  region  closer  to  the  WD  is  the 
boundary  layer  region  (discussed  further  in  §  10).  Standard  model 
temperatures,  for  the  M  values  we  are  considering,  are  so  high  at 
the  temperature  maxima  that  only  the  Rayleigh-Jeans  tail  contrib¬ 
utes  to  the  observed  spectra  from  those  regions.  (If,  in  addition,  as 
described  above,  half  of  the  total  liberated  potential  energy  were 
deposited  in  the  boundary  layer  region  the  Teff  there  would  be¬ 


Fig.  5. — Comparison  to  observed  spectra  of  model  with  17,000  K  accretion 
disk  face  and  associated  M  =  6.0  x  10“'^  Mgyr"*.  The  WD  mass  is  1.20  Mg.  The 
synthetic  (line)  spectmm  (heavy  line)  has  a  resolution  of  0.02  A.  The  synthetic 
spectrum  has  been  divided  by  4.5  x  1 0“*^.  The  corresponding  distance  to  QU  Car  is 
687  pc.  See  the  text  for  details. 

come  much  higher.)  The  bolometric  flux  contribution,  per  unit 
area,  of  the  innermost  annuli  far  exceeds  that  of  more  remote 
parts  of  the  accretion  disk.  Consequently  we  have  preserved  the 
standard  model  temperature  maxima  of  the  boundary  layer  re¬ 
gion  and  have  modified  the  temperatures  of  the  remaining  an¬ 
nuli.  This  choice  preserves  a  connection  with  a  particular  M 
through  the  sensitivity  of  the  total  accretion  disk  flux  to  the  con¬ 
tribution  of  the  highest  temperature  annuli.  For  convenience  of 
nomenclature,  we  designate  the  region  between  the  WD  and  the 
accretion  diskR=  1.3611RwDas  the  boundary  layer  even  though 
this  term  has  a  slightly  diflerent  connotation  in  the  literature. 

It  has  proved  possible  to  achieve  excellent  fits  by  adopting  an 
isothermal  profile  over  almost  all  of  the  accretion  disk  while  pre¬ 
serving  a  narrow  region  adjacent  to  the  boundary  layer,  discussed 
in  more  detail  below,  and  assigned  an  intermediate  temperature 
value.  The  narrow  region  is  important  in  achieving  a  fit  to  the 
short-wavelength  part  of  the  FUSE  spectrum.  The  highest  tem¬ 
perature  boundary  layer  region,  important  in  the  total  accretion 
disk  energy  budget,  contributes  to  the  flux  shortward  of  1000  A 
in  the  observed  spectrum  by  an  amount  that  varies  slowly  with  its 
temperature  (the  contribution  is  in  the  Rayleigh-Jeans  tail).  Chang¬ 
ing  the  boundary  layer  Teff  from  its  standard  model  value  produces 
appreciable  changes  in  the  system  flux  (thereby  disconnecting  from 
the  adopted  M  designation)  without  producing  a  sensitive  fit  to 
the  observed  FUV  flux.  The  narrow  intermediate-temperature  re¬ 
gion  provides  a  reasonably  sensitive  means  to  produce  a  good  fit 
to  the  observed  FUV  flux.  In  all  of  our  models,  the  contribution 
of  the  intermediate-temperature  region  near  1000  A  also  is  in  its 
Rayleigh-Jeans  tail. 

The  range  of  isothermal  values  that  we  have  tested  and  that 
produce  acceptable  fits  for  this  Mq  ,  is  limited.  The  range  extends 
from  15,000  K  at  M  =  1.0  x  10-’  Mg  yr^'  to  17,000  K  at  M  = 
1.0  X  10^®  Mq  yr^'  (higher  isothermal  temperatures  produce  a 
too-large  spectral  gradient). 

Figure  5  presents  a  sample  fit  to  the  observed  spectra,  again  for 
M  =  6.0  X  10  ^Mgyr  *.  The  isothermal  region  has  a  temperature 
of  1 7,000  K.  In  this  comparison  we  have  calculated  a  line  spectmm 
with  a  resolution  of  0.02  A.  The  resolution  of  the  FUSE  spectmm 
is  0.05  A.  Our  choice  of  synthetic  spectmm  resolution  was  to  permit 
an  accurate  correction  for  ISM  absorption,  described  in  §  8.  Figure  6 
shows  details  of  the  fit  to  the  FUSE  spectmm  and  the  FUV  end  of 
the  STIS  spectmm.  Extensive  absorption  lines  of  the  ISM  com¬ 
plicate  the  comparison.  Further  consideration  of  this  topic  is  in  §  8. 
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Fig.  6. — Detail  of  Fig.  5.  The  synthetic  spectrum  {heavy  line)  fits  the  peaks  of 
the  FUSE  spectrum  {gray  plot)  but  the  latter  spectrum  shows  strong  effects  of 
ISM  absorption.  The  fit  includes  the  short-wavelength  end  of  the  STIS  spectrum 
{light  line).  See  the  text  for  details. 

Figure  7  shows  the  fit  to  the  STIS  spectrum.  The  failure  to  fit 
Lya,  as  we  show  subsequently,  is  rectified  by  including  ISM  ab¬ 
sorption.  The  1200-1500  A  interval  is  alTected  by  the  ISM.  With 
the  exception  of  not  fitting  the  He  ii  emission  line,  the  fit  in  the 
1560-1710  A  interval,  shown  in  Figure  8,  generally  is  very  good. 
Finally,  Figure  9  shows  the  accurate  fit  to  the  red  end  of  the  SWP 
spectmm  and  the  blue  end  of  the  LWP  spectmm. 

The  bolometric  luminosity  for  the  M  —  6.0  x  10^^  Mq  yr^'" 
standard  model  accretion  disk  is  1.01  x  10^^  ergs  s^*.  The  bolo¬ 
metric  luminosity  of  the  empirical  model  is  3.95  x  10^®  ergs  s^*. 
Figures  5-9  are  representative  of  corresponding  plots  for  the  four 
largest  mass  transfer  rates  for  this  WD  and  secondary  star  mass. 

The  system  parameters  for  this  A/wq  are  in  Table  4.  The  default 
BINSYN  model  (before  imposing  a  nonstandard  accretion  disk 
temperature  profile)  assigns  standard  model  properties  (FKR92) 
at  the  various  annulus  radii,  including  the  annulus  semithickness. 
Those  properties  vary  slightly  with  M  but  the  tabulated  value  of  H, 
within  the  thickness  resolution  listed,  is  appropriate  for  the  range 
of  M  included  in  this  study.  Because  our  empirical  models  pre¬ 
serve  a  connection  to  corresponding  standard  model  M  values, 
the  H  value  remains  the  same  in  the  empirical  models. 

7.  QU  CARINAE  MODELS  WITH  Mwd  =  0.6  Mq,  q  =  0.83 

We  wish  to  test  the  possibility  of  acceptable  standard  model 
SED  fits  to  the  observed  spectra  over  the  full  range  of  possible 
Mwd  and  M  values.  Figure  2  indicates  that  the  smallest  possible 


Fig.  7. — Detail  of  Fig.  5  showing  fit  to  the  STIS  spectmm.  ISM  absorption 
affects  the  STIS  spectmm  shortward  of  1500  A.  The  lUE  contribution  has  been  sup¬ 
pressed  to  avoid  confusion  with  the  STIS  spectrum.  See  the  text  for  a  discussion. 


Fig.  8. — Expanded  scale  detail  of  Fig.  5,  synthetic  spectium  resolution  0.02  A. 
The  gray  plot  is  a  segment  of  the  STIS  spectmm,  resolution  0.1  A.  Note  that  the 
ordinate  scale  is  linear.  Except  for  a  few  details  and  the  He  ii  emission  line,  the  fit  is 
quite  accurate.  Note  the  doubling  of  the  He  ii  line.  The  synthetic  spectmm  normal¬ 
izing  factor  is  the  same  as  in  Fig.  5  (i.e.,  the  normalizing  factor  has  not  been  changed 
to  produce  a  better  fit  in  this  spectral  region). 

Mwd  for  this  system  is  about  0.6  Mq.  Again  from  Figure  2, 
A/wd  =  0.6  Mq  and  Mg  =  0.5  Mq  fall  on  the  i  =  55°  contour. 
From  Figure  4a  of  Panei  et  al.  (2000)  for  a  A/wd  =  0-6  Mq  with 
Teff  =  55,000  K,  we  find  Rwd  =  1-59  x  10^2  xhe  corre¬ 
sponding  radius  of  a  zero-temperature  WD  is  1.21  X  lO^^i?0.For 
models  with  a  range  of  M  values,  and  the  WD  =  55,000  K, 
the  tidal  cutoif  radius  of  the  accretion  disk  is  5  3R  wd  >  which  also  is 
0.840  Rq. 

We  tested  models  with  AT  =  0.3,  0.6,  1,  3,  6,  and  10  x 
10^^  Mq  yr^*  The  outer  part  of  a  standard  model  accretion  disk 
for  the  first  mass  transfer  rate  is  cooler  than  6000  K  and  is  unstable 
Smak(1982).  That  model  must  be  rejected.  In  all  remaining  cases, 
the  standard  model  produces  an  accretion  disk  profile  with  a  spec¬ 
tral  gradient  that  is  too  large.  The  effect  is  similar  to  the  ATwp  = 
1 .20  Mq  cases.  Figure  10,  illustrating  the  standard  model  for  AT  = 
3.0  X  10^^  Mq  yr^*,  is  representative  of  the  standard  model  syn¬ 
thetic  spectrum  fits.  As  with  the  AT  =  1 .20  Mq  case,  we  find  that 
no  standard  model  SED  fits  the  observed  spectra. 

Figure  11  presents  a  truncated  model  for  AT  =  3.0x 
lO^^AT0yr^'.  The  truncation  radius  isRtmnc  =  9.91Rwd-  The 
overall  fit  is  a  big  improvement  over  Figure  10  but  the  fit  to  the 
FUSE  spectmm  is  unsatisfactory.  The  same  type  defect  is  present 
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Fig.  9. — Detail  of  Fig.  5  showing  fit  of  the  synthetic  spectmm  {heavy  line)  to 
the  long-wavelength  end  of  the  SWP  lUE  spectmm  {light,  slowly  undulating  line) 
and  the  short- wavelength  end  of  the  LWP  lUE  spectmm  {light  jagged  line).  The  ac¬ 
curate  superposition  of  the  SWP  spectmm  on  the  red  end  of  the  STIS  spectmm 
{short,  heavy,  noisy  line)  is  visible  at  the  extreme  left  of  the  plot. 
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QU  Car  Model  System  Parameters,  Mwd  —  1-2  Mq 


Parameter 

Value 

Parameter 

Value 

^WD . 

1.20  MQ(adopted) 

ry(pole) . 

1.185  i?Q 

M, . 

1.2  Mq 

ry(point) . 

1.664  i?Q 

p . 

0.4542  days 

ry(side) . 

1.245  i?Q 

D . 

3.3275  «Q 

r^back) . 

1.348  i?Q 

flwD  . 

500.0 

logadpole) . 

4.39 

ft, . 

3.75 

logfldpoint) . 

-5.00 

4.30 

Uff.WD  . 

55,000  K  (nominal) 

log  3,  (back) . 

4.15 

rwD . 

0.00666  i?Q 

Fa . 

0.998  i?Q 

8.9 

0.00666 

H. . 

0.05  «o 

Note. — WD  refers  to  the  WD;  s  refers  to  the  secondary  star.  D  is  the  component 
separation  of  centers,  Qisa  Roche  potential,  specifies  the  outer  radius  of  the  ac¬ 
cretion  disk,  set  at  the  tidal  cutoff  radius,  and  is  the  accretion  disk  inner  radius. 
H  is  the  semiheight  of  the  accretion  disk  rim  (standard  model). 
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Fig.  11. — Fit  of  truncated  model  continuum  synthetic  spectrum,  WD  mass 
of  0.6  Mq,  to  combined  observed  spectra.  The  mass  transfer  rate  is  M  =  3.0  x 
10“^  Mq  yr“^  The  truncation  radius  isi^tnmc  =  9.91/?wd-  The  synthetic  spectrum 
fails  to  fit  the  observed  spectra  in  the  FUV.  Compare  with  Fig.  4. 


for  the  other  M  values.  We  conclude  that  no  truncated  model  pro¬ 
vides  a  SED  which  fits  the  observed  spectra  satisfactorily.  As  with 
the  Mwd  =  1.20Mq  case,  it  proves  possible  to  produce  good  qual¬ 
ity  synthetic  SED  fits  for  the  four  largest  mass  transfer  rates  we 
tested.  Figure  12  illustrates  the  result  for  M  =  3 .0  x  IQ-’Meyr-'. 
Fits  for  other  mass  transfer  rates  are  of  similar  quality.  Figure  13 
presents  the  FUSE  fit,  and  Figure  14  the  STIS  fit.  Figure  1 5  shows 
the  expanded  scale  fit,  1560-1710  A.  Figure  16  shows  the  fit  to 
the  lUE  spectra  and  the  accurate  connection  to  the  red  end  of  the 
STIS  spectrum.  Parameters  of  the  model  for  this  Mwd  ^re  in 
Table  5.  The  comments  describing  Table  4  also  apply  to  Table  5. 

In  summary,  over  the  full  range  of  permissible  Mwd  rnasses 
in  this  system,  no  acceptable  M  has  an  associated  standard  model 
SED  that  fits  the  observed  spectra  with  acceptable  accuracy.  The 
same  conclusion  follows  for  truncated  standard  models.  As  with 
the  Mwd  =1.2  Mq  case,  all  permissible  M  cases  have  an  asso¬ 
ciated  empirical  accretion  disk  model  that  provides  an  accurate 
fit  to  the  set  of  observed  spectra.  Further  discussion  is  in  §  10 

8.  INCLUSION  OF  A  MODEL  OF  THE  ISM 

We  have  developed  software  to  simulate  the  ISM  and  have 
used  this  routine  to  determine  ISM  properties  in  the  line  of  sight 
to  QU  Car.  A  custom  spectral  fitting  package  is  used  to  estimate 
the  temperature  and  density  of  the  interstellar  absorption  lines 
of  atomic  and  molecular  hydrogen.  The  ISM  model  assumes  that 


1000  1500  2000  _  2500  3000 
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Fig.  10. — Fit  of  standard  model  continuum  synthetic  spectrum  to  combined 
observed  spectra.  The  mass  transfer  rate  is  M  =  3.0x10  ^Mgyr  *,  andtheWD 
mass  is  0.6Mq.  The  synthetic  spectram  spectral  gradient  is  too  large.  Compare  with 
Fig.  3. 


the  temperature,  bulk  velocity,  and  turbulent  velocity  of  the  me¬ 
dium  are  the  same  for  all  atomic  and  molecular  species,  whereas 
the  densities  of  atomic  and  molecular  hydrogen,  and  the  ratios  of 
deuterium  to  hydrogen  and  metals  (including  helium)  to  hydro¬ 
gen  are  adjustable  parameters.  The  model  uses  atomic  data  of 
Morton  (2000, 2003)  and  molecular  data  of  Abgrall  et  al.  (2000). 
The  molecular  hydrogen  transmission  values  have  been  checked 
against  those  of  McCandliss  (2003).  Most  of  the  lines  in  the  FUSE 
spectrum  are  due  to  molecular  hydrogen. 

The  ISM  parameters  are  in  Table  6.  The  H  i  column  density  is 
in  reasonable  agreement  with  Drew  et  al.  (2003).  Although  Drew 
et  al.  present  evidence  for  a  two-component  ISM  toward  QU  Car, 
we  have  restricted  our  simulation  to  a  single  component.  The 
output  of  the  simulation  program  is  a  list  of  transmission  values 
at  the  same  spectral  resolution  as  the  system  synthetic  spectrum. 
The  product  of  the  ISM  transmission  function  and  our  0.02  A 
resolution  synthetic  spectrum  produces  a  corrected  synthetic 
spectrum.  We  have  used  the  M  =  6.0  x  10^’MQyr^^,MwD  = 
1 .20  Mq  case  to  illustrate  inclusion  of  the  ISM  simulation. 

The  ISM-corrected  synthetic  spectrum  line  density  is  so  great 
in  the  FUSE  region  that  it  is  not  meaningful  to  present  a  plot  com¬ 
parable  to  Figure  6  (and  a  color  plot  is  not  helpful).  Figure  17 
shows  a  partial  plot  of  the  FUSE  comparison;  note  the  cutoif  at 
912  A  as  compared  with  Figure  6.  The  ISM  lines  in  our  new 
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Fig.  12. — Fit  of  an  isothermal  model  synthetic  spectrum  {heavy  line)  to  com¬ 
bined  observed  spectra.  The  synthetic  spectrum  has  a  resolution  of  0.02  A.  The 
mass  transfer  rate  is  M  =  3.0  x  10  ^Mgyr  ^  and  the  WD  mass  is  0.6  Mq.  The 
synthetic  spectrum  has  been  divided  by  1 .57  x  10^^.  The  corresponding  distance 
to  QU  Car  is  406  pc.  See  the  text  for  details. 
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Fig.  13. — Detail  of  Fig.  12  showing  comparison  with  the  FUSE  spectrum 
{gray  plot)  and  the  short-wavelength  end  of  the  STIS  spectrum.  Note  the  strong 
ISM  absorption  lines  in  the  FUSE  spectrum. 

model  fit  the  observed  lines  well  in  many  cases,  including  line 
depths.  There  are  depth  discrepancies  in  some  cases,  and  our  model 
misses  some  observed  lines.  Our  ISM  model  is  a  reasonable  first 
order  representation  of  a  very  complex  ISM.  There  now  is  a  good 
fit  at  Lyy  and  Ly&  Large  unmodeled  residuals  remain  from  N  iv 
and  S  VI  (see  Fig.  1).  Hartley  et  al.  (2002)  discuss  features  like 
these  in  terms  of  a  wind  and/or  a  disk  chromosphere.  No  adjust¬ 
ment  of  parameters  in  our  present  set  of  models  can  produce  a  fit 
to  these  features,  and  they  represent  the  largest  remaining  resid¬ 
uals.  Note  the  differences  in  the  ordinates  in  Figure  1 7  and  Figure  1 . 
The  latter  apply  to  unreddened  fluxes  and  the  former  to  fluxes 
corrected  for  reddening. 

Figure  1 8  shows  the  fit  near  Ly/3.  Note  the  unmodeled  broad  O  vi 
features  at  1033  and  1037  A  and  the  improved  Ly/3  fit  as  compared 
with  Figure  6.  Figure  19  shows  the  fit  to  the  STIS  spectrum.  The 
fit  to  Lya  now  is  accurate  (compare  with  Fig.  7).  The  STIS  spec¬ 
trum  shows  a  number  of  unmodeled  absorption  features  includ¬ 
ing  the  N  V  line  pair  near  1 243  A  and  the  series  of  absorption  lines 
starting  near  1370  A. 

Although  we  have  illustrated  correction  for  the  ISM  for  a  sin¬ 
gle  empirical  model,  the  procedure  is  applicable  to  all  of  those 
models.  Our  ISM  model  fit  to  the  C  i  lines  in  the  FUSE  spectrum 
does  not  indicate  a  C  overabundance  in  the  ISM.  Drew  et  al.  (2003) 
present  evidence  of  a  C  overabundance  in  the  matter  undergoing 
transfer  from  the  secondary,  based  on  C  emission  lines. 
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Fig.  14. — Detail  of  Fig.  12  showing  comparison  with  the  STIS  spectrum. 
ISM  absorption  affects  the  STIS  spectrum  shortward  of  1500  A.  The  lUE  contri¬ 
bution  has  been  suppressed  to  avoid  confusion  with  the  STIS  spectrum. 
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Fig.  15. — Expanded  scale  detail  of  Fig.  12.  Note  that  the  ordinate  scale  is  linear. 
The  gray  line  is  a  segment  of  the  STIS  specffum,  resolution  0. 1  A.  The  synthetic 
spectrum  resolution  is  0.02  A.  Except  for  a  few  details  and  the  He  ii  emission  line, 
the  fit  is  quite  accurate.  Compare  with  Fig.  8.  The  synthetic  spectrum  normalizing 
factor  is  the  same  as  in  Fig.  12. 

9.  DETERMINATION  OF  SYSTEM  PARAMETERS 

Several  conclusions  follow  from  fhe  analysis  fo  this  point. 

1 .  No  standard  model  that  we  studied  fits  the  observed  spectra 
in  either  of  the  two  cases.  For  reasons  discussed  previously, 
we  believe  the  same  conclusion  would  be  tme  for  any  interme¬ 
diate  case,  and  we  believe  the  two  cases  are  limiting  cases. 

2.  No  truncated  standard  model  that  we  studied  produces  an 
acceptable  fit  to  observed  spectra  in  either  of  the  cases. 
Again,  for  reasons  already  discussed,  we  believe  that  situation 
would  be  true  for  any  intermediate  case,  and  that  the  two  cases 
are  limiting  cases. 

3.  Possible  M  values  are  bounded  for  small  values  by  an  un¬ 
stable  accretion  disk.  This  limit  is  about  M  =  l.Ox  10^^  A/q  yr^*. 

4.  A  possible  constraint  at  high  M  is  theoretical:  expansion  of 
the  WD  to  a  red  giant  and  consequent  engulfing  of  the  entire  bi¬ 
nary  (Paczyhski  &  Zytkow  1978;  Sion  et  al.  1979).  However,  the 
studies  cited  considered  spherical  accretion  and  do  not  apply  to 
disk  accretion  (Paczyhski  &  Zytkow  1978).  The  effects  of  mass 
transfer  onto  WDs  via  accretion  disks  have  been  studied  by  Piro 
&  Bildsten  (2004)  based  on  a  spreading  layer  model.  These  au¬ 
thors  consider  M  values  up  to  1.5  x  10^^  Mq  yr^*  and  find  no 
significant  expansion  of  the  WD.  Note  that  their  study  considers 
the  mass  actually  transferred  to  the  WD,  and  in  the  case  of  QU  Car 
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Fig.  16. — Detail  of  Fig.  12  showing  the  fit  of  the  synthetic  spectrum  {heavy 
line)  to  the  long-wavelength  end  of  the  SWP  lUE  spectrum  {light,  slowly  un¬ 
dulating  line)  and  the  short-wavelength  end  of  the  LWP  lUE  spectrum  {light  jagged 
line).  The  accurate  superposition  of  the  SWP  spectrum  on  the  red  end  of  the  STIS 
spectrum  {short,  heavy,  noisy  line)  is  visible  at  the  exffeme  left  of  the  plot. 
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TABLE  5 

QU  Car  Model  System  Parameters,  Mwd  =  0-60  Mq 


Parameter 

Value 

Parameter 

Value 

^WD . 

0.60  MQ(adopted) 

ry(pole) . 

0.8748  itg 

M, . 

0.50  Mg 

ry(point) . 

1.235  Rq 

p . 

0.4542  days 

ry(side) . 

0.9169  itg 

D . 

2.56556  ttg 

r^back) . 

0.9977  itg 

t^WD  . 

162.2 

log&(pole) . 

4.27 

n, . 

3.47364 

log  appoint) . 

0.43 

55°  (adopted) 

4.19 

Tjff.WD  . 

55,000  K  (nominal) 

log&(back) . 

4.03 

rwD . 

0.0159  Rq 

Fa . 

0.84  itg 

logSwD  . 

7.8 

rb . 

0.0159  itg 

H . 

0.05  itg 

Notes. — WD  refers  to  the  WD;  s  refers  to  the  secondary  star;  D  is  the 
component  separation  of  centers,  fi  is  a  Roche  potential;  specifies  the  outer 
radius  of  the  accretion  disk,  set  at  the  tidal  cutoff  radius;  and  is  the  accretion 
disk  inner  radius.  H  is  the  semiheight  of  the  accretion  disk  rim  (standard  model). 

there  is  a  wind  which  may  carry  away  an  appreciable  mass 
(Hartley  et  al.  2002).  At  =  1.0  x  10^®  Mq  yr^'  and  larger 
Piro  &  Bildsten  (2004)  point  out  that  hydrogen  is  expected  to 
be  burning  steadily  on  the  WD  surface,  leading  to  luminosities 
in  the  range  lO^^-lO^^  erg  s^^  A  luminosity  that  large  exceeds 
the  disk  luminosity  and  so  would  be  incompatible  with  the 
observed  double  14686  emission  line.  We  take  M=  1.0  x 
10^®  Mq  yr^^  to  be  an  upper  limit  to  acceptable  mass  transfer 
rates. 

5.  Nonstandard  models  can  be  found  that  produce  acceptable 
fits  to  the  combined  observed  spectra.  We  achieved  successful  fits 
with  a  hot  accretion  disk  region  bracketed  between  the  WD  and  a 
radiusR  =  1.361  IRwd  (which  we  designate  as  aboundary  layer), 
an  intermediate-temperature  region  adjacent  to  the  hot  region,  and 
an  isothermal  region  extending  to  the  tidal  cutoff  boundary.  The 
hot  region  Teff  was  set  equal  to  the  standard  model  temperature 
maximum  to  preserve  a  designation  connection  with  the  standard 
model  M.  Because  the  peak  emissions  of  both  the  hot  region  and 
the  intermediate-temperature  region  are  at  much  shorter  wave¬ 
lengths  than  the  short-wavelength  limit  of  the  FUSE  spectmm,  we 
have  no  good  way  to  impose  tight  temperature  constraints  on  those 
regions. 

Table  7  compares  properties  of  models  that  fit  the  observed 
spectra  well.  Columns  (3),  (4) ,  and  (5)  list  the  values  for  the 
boundary  layer,  the  intermediate  region,  and  the  isothermal  re¬ 
gion.  The  last  four  lines  list  the  possible  range  of  isothermal  re¬ 
gion  Teff  values.  The  first  four  lines,  for  Mwd  =  1  -2  Mq,  present, 
first,  a  comparison  of  line  1  with  line  6  for  the  same  Tefif  values 
of  the  intermediate  and  isothermal  regions.  Line  6,  for  Mwd  = 
0.6  Mq,  corresponds  to  a  slightly  smaller  accretion  disk  (com¬ 
pare  in  Table  4  and  Table  5)  and  the  value  of  i  for  Table  4  is  40° 
versus  55°  for  Table  5.  Line  3  and  line  4  preserve  the  same  value 


TABLE  6 

Parameters  of  ISM  model 


Parameter 

Value 

A'tH) . 

2.5  X  10^"  cm-2 

D/H . 

1  X  10-5 

A'(H2)  . 

7  X  10‘5  cm-2 

Velocity . 

0.0  km  s“' 

Temperature . 

350  K 

Metallicity . 

0.01  (1.0  =  solar) 

910  920  930  940  950  960  970  980 

Wavelength  (A) 


Fig.  17. — ISM-corrected  version  of  Fig.  6.  The  spectral  resolution  of  the 
synthetic  spectrum  is  0.02  A.  Strong  absorption  lines  of  N  iv  and  S  vi  {gray  under¬ 
lying  plot)  remain  unmodeled,  possibly  from  a  wind  or  accretion  disk  chromo¬ 
sphere.  See  Fig.  1  for  identifications.  There  now  is  a  good  fit  to  Ly7  and  Ly(5.  Note 
the  linear  ordinate  scale.  See  the  text  for  a  discussion. 

of  isothermal  region  Tejf  and  illustrate  the  increase  of  bolometric 
luminosity  with,  primarily,  the  change  in  the  boundary  layer  Teff. 
Note  that  the  calculated  distance  to  QU  Car  is  the  same  for  these 
two  lines,  indicating  that  the  normalizing  factor  to  fit  the  synthetic 
spectrum  to  the  observed  spectra  depends  strongly  on  the  isother¬ 
mal  region  and  is  fairly  insensitive  to  the  boundary  layer  Tefif- 
Column  (6)  lists  the  bolometric  luminosity  of  the  standard  model 
for  comparison  with  the  bolometric  luminosity  of  the  empirical 
model  in  column  (7).  Note  that  the  entries  in  the  two  columns  are 
roughly  comparable,  in  agreement  with  oiu  earlier  discussion  on 
preserving  a  cormection  between  the  empirical  model  and  the  stan¬ 
dard  model  M.  The  last  column  lists  the  distance  to  QU  Car  for 
each  empirical  model.  The  range  of  distances  is  not  large  for 
Mwd  =  1.2  Mq,  and  all  are  in  rough  accord  with  the  Hipparcos 
distance  (§  2).  The  smallest  distance  for  Mwd  =  0.6  Mq  (332  pc) 
is  in  marginal  disagreement  with  the  Hipparcos  distance. 

There  is  no  clear  distinction,  among  the  acceptable  models,  of 
the  quality  of  fit.  It  is  important  to  note  that  the  FUSE,  STIS,  and 
lUE  spectra  were  obtained  at  different  times,  that  they  were  fitted 
together  by  visual  estimate,  and  that  the  QU  Car  luminosity  varies 
on  a  short  timescale  (Schild  1969;  Knigge  et  al.  1994).  Because  of 
this  continuing  variation,  a  new  study  with  new  data  likely  would 
find  larger  differences  tfom  the  present  data  than  the  residuals  in 
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Fig.  1 8. — Comparison  of  ISM-corrected  synthetic  spectrum  and  FUSE  spec¬ 
tmm  near  Ly/?.  The  ISM  correction  provides  a  good  fit  to  a  number  of  the  narrow 
ISM  lines.  Note  the  linear  ordinate  scale.  Note  the  remaining  unrepresented  broad 
O  VI  lines  at  1033  and  1037  A.  These  features  probably  associate  with  a  wind  or 
chromosphere  not  included  in  the  model. 
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Fig.  19. — Synthetic  spectrum  fit,  including  an  ISM  correction,  to  STIS  spec¬ 
trum.  The  heavy  line  is  the  synthetic  spectrum,  and  the  gray  line  is  the  observed 
specttnm.  The  lUE  contribution  has  been  suppressed  to  avoid  confusion  with  the 
STIS  spectrum.  There  now  is  a  good  fit  to  Lya,  in  contrast  to  Fig.  7,  but  unmodeled 
features  remain.  See  the  text  for  a  discussion.  Note  the  linear  ordinate  scale. 

our  present  model  fits.  Although  accurate  fits  in  some  spectral  re¬ 
gions  are  possible  (e.g.,  Fig.  8),  the  fit  accuracy  should  not  be 
overinterpreted. 

The  fit  in  the  920-1200  A  region  is  fairly  sensitive  to  the  ^eff 
of  the  intermediate-temperature  region.  The  width  of  that  region 
arbitrarily  was  set  to  two  (BINSYN)  annuli,  and  their  corre¬ 
sponding  widths  depend  on  the  choice  of  45  annulus  divisions  of 
the  accretion  disk  in  the  BINSYN  model.  By  eye  estimate,  the 
fits  with  the  present  models  are  optimal.  We  suspect  that  a  choice 
of  one  of  the  models  over  another,  based  on,  say,  a  test,  could 
arise  from  numerical  eflfects  (interpolation  approximations,  etc.) 
rather  than  from  a  closer  representation  of  physical  reality.  Since 
the  largest  remaining  residuals  arise  from  unmodeled  wind  or 
chromospheric  features,  these  models  present  the  closest  repre¬ 
sentation  of  the  QU  Car  data  currently  available. 

Tables  4  and  5  list  the  orbital  parameters  used  for  the  two  val¬ 
ues  of  Mwd-  Since  there  is  no  light  curve  to  simulate,  the  of 
the  secondary  star  is  of  importance  only  to  the  extent  that  it  could 
affect  the  long-wavelength  end  of  the  system  synthetic  spectrum. 
A  nominal  6500  K  secondary  in  the  A/wd  =1.2  Mq  system  pro¬ 
vides  no  evidence  of  contamination  in  the  system  synthetic  spec¬ 
trum.  Initial  retention  of  that  Tejf  in  the  Mwd  =  0.6  Mq  system, 
followed  by  substitution  of  a  3500  K  secondary  produces  no  de¬ 
tectable  effect  on  the  system  synthetic  spectrum. 

We  are  confident  that  Mwd  masses  could  be  selected  between 
our  limiting  values  of  0.6  and  1 .2  Mq,  combined  with  appropriate 


values  of  M^  and  i  from  Figure  2,  assigned  a  range  of  M  values  and 
produce  SED  fits  to  the  observed  spectra  comparable  in  quality  to 
Figure  5  or  Figure  12.  Table  7  represents  the  extremes  of  possible 
system  models,  based  on  existing  spectroscopic  data.  Because  of 
its  consistency  with  the  empirical  mass-period  relation  (Warner 
1 995,  eq.  [2. 1 00]),  we  favor  the  Afwo  =1-2  Mq  solution. 

10.  DISCUSSION 

We  have  appropriated  the  term  “boundary  layer”  to  desig¬ 
nate  the  region  between  the  WD  and  the  accretion  disk  R  — 
1 .361  li?wD-  This  region  more  generally  is  associated  with  a  very 
hot  region  occupying  the  same  space  but  predicted  to  emit  hard 
X-rays  if  optically  thin  or  soft  X-rays  if  optically  thick  (Warner 
1995,  §  2.5.4).  This  prediction  is  modified  if  a  wind  carries  off 
part  of  the  accretion  mass  (Hoare  &  Drew  1991;  Hoare  &  Drew 
1993).  In  the  latter  case,  predicted  boundary  layer  values  of  Teff 
are  no  higher  than  the  values  adopted  in  our  models.  Based  both 
on  its  geometric  location  and  the  much  higher  local  temperature 
than  in  the  rest  of  the  accretion  disk,  our  appropriated  term  ac¬ 
tually  is  in  keeping  with  the  usual  sense  of  the  term. 

Adoption  of  the  intermediate-temperature  layer  obviously  is  a 
cmde  approximation.  A  temperature  discontinuity  at  its  inner  and 
outer  boundary,  as  in  our  present  models,  is  not  physically  cred¬ 
ible.  Flowever,  this  model  is  a  much  better  fit  to  observation  than 
with  our  earlier  experimental  (§  6)  models,  which  adopted  a  sharp, 
nonetheless  more  gradual  temperature  drop  to  the  isothermal  disk 
region  than  in  this  paper. 

The  outermost  accretion  disk  region  contributes  little  to  the 
system  synthetic  spectrum.  The  Tefif  values,  for  a  limited  region 
which  would  vary  from  one  M  to  another,  could  be  closer  to  stan¬ 
dard  model  values  than  the  isothermal  assumption  without  pro¬ 
ducing  a  detectable  change  in  the  system  synthetic  spectrum.  On 
the  other  hand,  it  is  known  that  temperatures  in  the  outer  region 
of  an  accretion  disk  need  an  upward  correction  from  the  standard 
model  to  account  for  impact  heating  (Lasota  2001)  and  tidal  heat¬ 
ing  (Buat-Menard  et  al.  2001). 

The  accretion  disk  temperature  profiles  found  in  this  study  dif¬ 
fer  from  profiles  found  for  other  NL  systems.  In  both  SDSS  J0809 
(Linnell  et  al.  2007a)  and  IX  Vel  (Linnell  et  al.  2007b)  the  final 
model  included  an  iimer  isothermal  region  with  an  outer  region  that 
followed  a  standard  model.  In  the  case  of  MV  Lyr  (Linnell  et  al. 
2005)  in  a  high  state,  the  final  model  had  a  standard  model  accretion 
disk  extending  from  an  inner  truncation  radius  to  an  intermediate 
radius  with  an  isothermal  radius  beyond.  It  is  of  interest  that  all 
of  the  NL  systems  included  in  these  studies  have  accretion  disks 
whose  temperature  profiles  clearly  differ  from  the  standard  model. 


TABLE  7 

Properties,  Empirical  Mass  Transfer  Models 


M 

(Mo) 

(1) 

M 

(Mo  yr-') 

(2) 

7;ff(K) 

Aol  (ergs 

s-‘) 

D 

(pc) 

(8) 

Boundary  Layer 

(3) 

Intermediate  Region 

(4) 

Isothermal  Region 

(5) 

Standard  Model 
(6) 

Empirical  Model 

(7) 

1.20.. 

1.0  X  10“^ 

185,000 

90,000 

15,000 

1.68  X  10^'^ 

1.06  X  10” 

550 

1.20.. 

3.0  X  10“^ 

240,000 

100,000 

16,000 

5.05  X  10^'^ 

2.18  X  10” 

610 

1.20.. 

6.0  X  10“^ 

290,000 

110,000 

17,000 

1.01  X  10” 

3.95  X  10” 

687 

1.20.. 

1.0  X  10“'^ 

330,000 

110,000 

17,000 

1.68  X  10” 

5.52  X  10” 

687 

0.60.. 

1.0  X  10“^ 

100,355 

70,000 

14,000 

2.56  X  10” 

2.77  X  10” 

332 

0.60.. 

3.0  X  10“^ 

132,000 

90,000 

15,000 

7.68  X  10” 

7.11  X  10” 
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11.  SUMMARY 

The  analysis  in  this  study  leads  to  the  following  conclusions: 

1 .  Standard  model  SEDs,  within  the  tested  ranges  of  M  and 
Mwd>  are  too  hot  to  fit  the  combined  set  oiFUSE,  STIS,  and  lUE 
spectra  of  QU  Car  with  acceptable  accuracy. 

2.  Truncated  standard  model  SEDs  fail  to  fit  the  combined  set 
of  FUSE,  STIS,  and  lUE  spectra  of  QU  Car  for  similar  reasons. 

3 .  Alternative  nonstandard  models  have  SEDs  that  accurately 
fit  the  set  of  observed  spectra.  These  models  individually  have  a 
boundary  layer  equal  to  the  maximum  accretion  disk  of 
an  associated  standard  model  M.  A  narrow  adjustable  Tefp  annu¬ 
lus  is  needed  adjacent  to  the  boundary  layer  to  “fine  tune”  the 
SED  fit  to  the  FUSE  spectrum.  The  remainder  of  the  accretion 
disk  is  isothermal  to  within  the  accuracy  of  the  SED  fit.  The  model- 
to-model  residuals  are  so  similar  that  it  is  not  possible  to  choose  a 
best  model.  The  QU  Car  WD  mass  is  greater  than  0.6  Mq  and  less 
than  the  Chandrasekhar  limit;  a  solution  with  Mwd  =12  Mq  is 
consistent  with  an  empirical  mass-period  relation.  The  M  from  the 
secondary  star  is  between  1.0  x  10^^  and  1.0  x  10^^  A/0yr^^  The 
set  of  models  agree  on  a  limited  range  of  possible  isothermal  region 
Teg-  values  between  14,000  and  18,000  K. 

The  orbital  inclination  is  between  i  =  40°  and  i  =  60°.  The 
Hipparcos  distance,  610  pc,  is  representative  of  the  empirical 
model  results. 

4.  A  model  of  the  ISM  combines  with  the  empirical  accretion 
disk  models  to  provide  an  improved  fit  to  the  observed  FUSE 
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and  STIS  spectra.  The  largest  remaining  residuals  are  from  un¬ 
modeled,  broad,  high-excitation  features  probably  associated  with 
a  wind  or  a  disk  chromosphere. 

5.  The  accretion  disk  dominates  the  system  flux  in  the  FUSE, 
STIS  ,  and/L'E'  spectral  regions.  An  assumed  55,000  K  WD  pro¬ 
vides  a  negligible  flux  contribution,  as  does  the  accretion  disk  rim 
and  the  secondary  star. 
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